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This paper presents the theory of directional couplers using matrix 
formulation for lumped, distributed, or any linear time-invariant black-box. 
The starting point is the matrix theory of uniform multiple coupled trans- 
mission lines in terms of which are brought in the concepts of characteristic 
impedance matrix, propagation matrix, reflection matrix, and scattering 
matrix. Next, we use the results obtained with the theory of multiple coupled 
transmission lines to derive expressions for the loading impedance and 
voltage ratios for distributed directional couplers. We do this using the 
spectral properties of the matrices. Then we generalize the concepts of 
characteristic impedance matrix and propagation matrix for a class of 
black-boxes with the aid of the A, B, C, D transmission matrix. We give 
conditions for the loading impedance and expressions for the voltage 
ratios, using the spectral theory of the transmission matrix. We discuss 
the physical significance of the directional coupler effect at all frequencies 
in a vector-matrix framework and analyse in detail some lumped direc- 
tional couplers. Finally we discuss hybrid (lumped and distributed) 
directional couplers. 

I. INTRODUCTION 

The directional coupler is an important device in many transmission 
systems. The theory for the electrical design of certain types of dis- 
tributed-parameter directional couplers is well established through the 
contributions of a number of researchers. 1-7 

The purpose of this paper is to present the theory of transmission-line 
symmetric directional couplers in matrix form and then extend that 
theory to arbitrary lossless reciprocal circuits. There are many nota- 
tional and conceptual advantages in using a matrix formulation which 
give further insight for the synthesis of different types of directional 
couplers. 
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The starting point is to examine the coupled-line directional coupler 
as a particular application of multiple coupled transmission line theory 
and then to use the concepts of characteristic impedance matrix and 
propagation matrix to examine the whole problem, rather than one 
line at a time or one mode at a time. By using a matrix approach the 
fundamental properties of the modes become more evident and thus 
physical intuition and mathematical reasoning blend to give a clearer 
picture of the situation. For instance, when the lines are no longer 
identical it is the eigenvectors of the matrices which provide the in- 
formation of how to extend the mode concept. 

II. MULTIPLE COUPLED TRANSMISSION 



2.1 Propagation and Characteristic Impedance Matrices 

Several researchers 8 " 10 have analyzed the behavior of a set of multiple 
coupled transmission lines. For reference in subsequent sections of this 
paper we present a brief account of this theory. For simplicity the 
discussion is restricted to two identical lines operating in the TEM 
mode. Fig. 1 shows schematically a differential section of two lines 
and a ground plane. The circuit of which Fig. 1 is a differential section 
obeys the following vector differential equations* in the steady state 
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By solving for I in (1) and substituting its value in (2), or solving 
for V in (2) and substituting its value in (1) the following are obtained: 



d 2 V 

V, = ZYV, 

dx 



(3) 



* The matrices V, I, Z, Y are all functions of frequency. For simplicity in the 
notation this dependency is not explicitly indicated. 
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Fig. 1 — Differential section of two coupled lines above a ground plane. 
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Define the matrices r, Y and Z as follows: 

r 4 VZY, 
Y = z'r, 
Z = Y fJ . 
In terms of the matrices defined, the solution of (3) and (4) is* 
V(,-) = e" r *v + + e ri v_ , 
l(x) = Yo(e" r 'v + - e rx v_), 



(4) 

(5) 

(6) 
(7) 

(8) 
(9) 



where the 2-vectors v+ and v_ are arbitrary constants (dependent on 
frequency) which depend on the boundary conditions. Their inter- 
pretation is very similar to the one of single line theory, v+ may be 
called the forward wave and v_ the reflected wave. 

Equations (3), (8) and (9) have the same form as the single trans- 
mission line equations. For this reason the matrices r, Y , Z , are 
called propagation matrix, characteristic admittance matrix, and char- 
acteristic impedance matrix, respectively. Many properties of the 
single transmission line hold for the multiple case. In particular, if 
the set of lines is terminated in a network whose open circuit impedance 
matrix is equal to the characteristic impedance matrix of the set of 
lines, a vector of incident voltage waves traveling down the lines will 
experience no reflection. The manner in which the voltages and currents 



* In evaluating VZY to calculate r the convention is made to associate with r 
eigenvalues whose real part are positive and with — r the ones with negative real 
part. 
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in the lines interact is best understood by examining the matrix func- 
tions e~ r ' and e Ts . (See Appendix.) 

2.2 Reflection and Scattering Matrices 

In single transmission line theory the voltage reflection coefficient 
T R at a discontinuity is defined as 

v- = T R v + , (10) 

where v + is the incident voltage and v. is the reflected voltage. The 
expression for T R in terms of the characteristic admittance Y of the 
line and the input impedance Z L at the discontinuity is 



r« = 



Zr,Y n ~ 1 



(ID 



Z L Y + 1 

The concept of reflection coefficient can be generalized very simply 
for the case of a multiple set of lines. Consider Fig. 2 depicting a pair 
of coupled fines of length I and matrices r, and Y terminated at 
x = in a device of open circuit impedance matrix Z L . According 
to (8) and (9) for x = 

V(0) = v + + v_ , (12) 

1(0) = Y (v + - v_). (13) 

The box marked Z L obeys 

V(0) = Z L I(0). (14) 

Substitution of V(0) and 1(0) from (12) and (13) into (14) gives 

v + + v_ = Z*Y (v + - v_), (14b) 

from which solving for v_ 

v_ = (Z L Y + J)-f(Z L Y 9 - I)v + . (15) 

I r,Y i 



i 

X = 0! 







Fig. 2 — Pair of coupled lines terminated in a circuit of impedance matrix Zl. 
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Thus, defining the reflection matrix T R as a generalization of (10) 
according to 

v_ = r R v + , (16) 

it follows from (15) that 

T R = (Z L Y + IT'CZtY. - I), (17) 

which is the generalization to multiple lines of (11). 

The scattering matrix S is a reflection matrix in which the incident 
and reflected voltages are normalized.* It is defined by 

t- = Sv + ; (18) 

where 

v + = Y»v + , v_ = Y*v_ , (19) 

V = YjV, I = Zft. (20) 

From (14b) it follows that 

Z*(v + + v_) = Z,Y [Z»(v + - y-)]. (21) 

Defining the normahzed load impedance matrix Z L according to 

Z L = Y*Z,Y* , (22) 

equation (21) may be rearranged as follows: 

v_ = (Z L + iy\Z L - I)v + . (23) 

Comparison of (23) and (18) gives 

S = (Zl + 1)-\Zl - I). (24) 

From either (17) or (24) it is clear that if the load has an open circuit 

impedance matrix equal to Z of the lines the reflected wave is zero 

since both the reflection matrix r R and the scattering matrix S vanish. 

The complex power in terms of the normahzed variables is given by 

p = y + V = M. (25) 

(The superscript -f denotes transposed conjugate of a matrix.) For a 
lossless device the incident power must be equal to the reflected power, 



* In the literature on the scattering matrix the variables are normahzed with 
respect to a diagonal matrix (usually a real matrix). Here the matrix may be complex 
ana not necessarily diagonal. This derivation provides the physical interpretation 
for a scattering matrix normalized with respect to a complex nondiagonal matrix. 
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hence, from (18) 

p = r.v. = v :s + s v+ = r + * + , (26) 

which means that 

S^S - I, (27) 

that is, the scattering matrix is unitary. 

2.3 Terminal Matrices of a Set of Lines 

Viewed from its terminal behavior the two coupled lines and ground 
of Fig. 3 can be studied as a four port which may be characterized by, 
among others, an impedance matrix C or a transmission E matrix. 
(The E matrix is the extension to 2N ports of the concept of the A, B, 
C, D parameter matrix of two ports. See Ref. 11.) 

The < matrix may be written in partitioned form as follows : 

(csch[r-qZ )' ~| 
coth (r-Z)Z J' 



C = 



f coth (r-QZg 
[_csch(r-Z)Zo 



(28) 



(The prime indicates the transpose matrix.) The E matrix, written in 
partitioned form, is 

cosh r • I | sinh ( r-l)Z ~\ 

r'J J" 



F __ ["a I bi _ r_cps 

l_c I dJ " Ly si 



sinh r-l\ cosh r 



(29) 



If the two lines in Fig. 3 are identical then r is symmetric and hence 
r' = r so that in (29), A = D. Furthermore, for this case all the 
matrices Z , Y„ , r, A, B, C, D and their analytic functions commute 
and therefore may be treated unambiguously as ordinary numbers. 
For instance (29) may be written 



E = 



cosh r • I 


Zosinh r-l 


-=-svah r-l 
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Fig. 3 — Two coupled lines above a ground plane forming a four-port. 
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The analogy between the E matrix of (30) and the A, B, C, D parameter 
matrix of a single transmission line considered as a two port is now 
complete. 

III. SYMMETRICAL TRANSMISSION DIRECTIONAL COUPLERS 



3.1 Derivation of Matching Impedances 

A symmetric matched directional coupler is a four port device whose 
scattering matrix is of the form of S of (31) when the ports are ad- 
equately numbered. /3 and 5 are, in general, frequency dependent. 
It is well known in microwave circuits that, given any lossless reciprocal 
4-port, if all ports are matched, then the device is a directional coupler. 12 
By properly numbering the ports it may be assumed that the directional 
coupler has the following scattering matrix 



S - 
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(31) 



For the case of two identical lossless coupled lines equally loaded 
at the four ports all that is necessary to obtain a directional coupler 
is to match one of the ports. This results from the great symmetry. 

Fig. 4 shows two identical lines loaded at ports 2, 3, and 4 with 
equal resistances, R. To calculate the driving point impedance at 
port 1 the following procedure can be used. Consider Fig. 5. The 
equation 



Z = 



flA+ B 
RC + D 



Z 2 i ^22J 



(32) 




7777777777777777777777777777777T 
2 4 



Fig. 4 — Two-terminal circuit formed by a pair of coupled lines above a ground 
plane loaded with resistances at three ports. 
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is the matrix version of the well known formula for calculating the 
driving point impedance at one port of a two-port when the other 
port is loaded. In (32) Z is the open circuit impedance matrix of the 
two-port at the left in Fig. 5, A, B, C, D are the matrices given by (29). 
Once the Z of (32) is known the resistor R may be connected to 
port 2 as shown in Fig. 6 and the driving point impedance at port 1 
calculated by a second application of (32) although now, instead of 
matrices, A, B,C, D are scalars. The matrices in (32) are all of the form 



K = 



a. 8 



(33) 



hence the methods of the appendix are applicable. Using (30) and (32) 



Z = R,p + R 2 g 



__ [ to + g) 
" U(P - q) 
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where Ri and R 2 form the spectral set of Z and are 
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and p and g are the eigenvalues of Z and are 

cosh y + 1 + -=p- sinh y + 1 

V =Rr » 

-S= sinh 7 + 1 + cosh 7 + 1 

Z~ 

cosh7 - Z + ^-sinh7~Z 

!-»j " 

■=r sum 7 I + cosh y I 

The symbols y + , y~ denote the eigenvalues of T; and Z* , Z~ the 
eigenvalues of Z . 
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Fig. 5 — Intermediate circuit for calculating the driving point impedance at 
port 1 of the circuit in Fig. 4. 
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Fig. 6. — Second intermediate circuit for calculating the driving point im- 
pedance at port 1 of the circuit in Fig. 4. 



To transform the Z,, of (32) to A, B, C, D matrices the well known 
formulas : 



A — Z n Z 2 i , 

15 = ld\\L 21^22 "12 ) 

C = Z21 , 
D = Z 21^22 j 

which hold for matrices (and hence for scalars considering them as 
1 X 1 matrices) may be used. The second application of (32) with 
the aid of (38)-(41) gives 



(38) 
(39) 
(40) 
(41) 



_ (p + q)R + 2pq 
2R + v + ? 



(42) 



For port 1 to be matched z must equal R, hence, the condition for 
directional coupler effect is 



which reduces to 



_ (p + q)R + 2pq 
2R + p + q ' 



R = Vpq. 



(43) 



When the values given by (36) and (37) are substituted in (43), after 
some algebra, the following equation results : 



SF " z§^) sinh 7+z sinh y ~ l + (§■ ~ tJ sinb y ~ l cosh T+z 



+ (ff- - ff) smh T + Z cosh 7 ~Z = 0. 



(44) 



If a matched directional coupler at all frequencies is desired, (44) 
must be satisfied at all frequencies. Some possible mathematical solu- 
tions are the following: 
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(i) Make the three terms in parentheses in (44) vanish. This happens if 

Z* = Z~ = R. (45) 

(ii) A second possibility is to make 

7 + = 7~ = 7- (46) 

Then (44) reduces to 

+ (f" + f" " h " §=) sinh yl cosh yl = °- (47) 

To make the quantities in the parentheses in (47) vanish, R is given by 

R = VZ + Zo . (48) 

Equation (45) or (46) and (48) give possible conditions for directional 
coupler effect at all frequencies. If a narrow band directional coupler 
is desired one may match the coupler at the discrete frequencies which 
satisfy (47). Since (47) is a transcendental equation it is not unreason- 
able to expect an infinity of roots. For instance, suppose (46) holds, 
and I is made so that 

sinh 7^ = 0. (49) 

Then the device will be a directional coupler at the frequencies that 
are roots of (49). 

Some explicit relationships for two identical lossless lines of inductance 
per unit length L xl , mutual inductance per unit length L 12 , capacitance 
per unit length of one line alone C, and capacitance between the two 
lines per unit length C M , are : 

7 + = j V(L„ + L 12 )C o>, (50) 



7" = jViL^ -L 12 )(C + 2C^)co, (51) 

Zt = / (Lu + Llz) , (52) 

V(L n + L 12 )C 

ZZ = ten ~ L ") (53) 

V(L„ - L 12 )(C + 2C M ) 
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The condition expressed by (45) implies, equating (52) and (53) 

L u C + C M 



l-l\2 Ci 



(54) 



C + Cm 

While the condition expressed by (46) implies 



R = ^P^ 



and (48) reads 


L u C + C„ 

■^12 C M 


R 1 L » 


R \c + c M 



(56) 



which is the same as (55). Equation (54) implies that the lines have 
negative mutual inductance, which is not achievable with parallel 
lines. This condition can, however, be satisfield with counter-wound 
lumped elements. 

If (56) holds, which implies (46) then (49) implies 

ft -0,1, 2, .... (57) 



lVL 12 C M - L n (C + CZ) ' 

That is, for the frequencies given by (57), independent of the value 
of the loads (as long as they are all equal) the lines will be matched 
and will exhibit directional coupler effect. 

3.2 Frequency Dependency of the Coupling Between the Ports 

Equations (54) and (55) or alternatively (56) and (55) are not fre- 
quency-dependent. This means that the resulting circuit will be matched 
for all frequencies. Therefore, the directional coupler effect will exist 
for all frequencies, meaning that the coupling between uncoupled ports 
is zero at all frequencies. However, the coupling between coupled 
ports is frequency-dependent. This dependency is derived as follows. 

Considering the coupled transmission lines as the load to four un- 
coupled lines, each of characteristic impedance R as shown in Fig. 7, 
the scattering matrix of the load is calculated according to (24) and 
(22) with 



Y„ - 
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1 
R. 



(58) 



and Z L given by < of (28). Applying the methods of the appendix the 
matrix £ may be expressed as follows: 



C = RaXfa + RbXfC + RcXfc + RdXfd , 



(59) 



where R„ , R 6 , R„ , R d are the spectral set of C and they are given by 
(165)-(168) of the appendix and X fo , \ {h , X fe , X f(J are the eigenvalues 
of C and are given by 



X ro = coth (y + l)Z + Q + cosch (y*l)Z* , 
X f6 = coth (y + l)Z + - csch (y*l)Z+ , 
X fc = coth (y~l)Z~ + csch (y"l)Z~ ? , 
X fe = coth (y~l)Z~ — csch (y~l)Z" . 



(60) 
(61) 
(62) 
(63) 
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COUPLED LINES OF 
IMPEDANCE MATRIX 

Z L =? 



Fig. 7 — Two coupled lines above a ground plane of impedance matrix given 
by equation (28) serve as load to four uncoupled lines of characteristic imped- 
ance R for calculation of scattering matrix of coupled lines. 



DIRECTIONAL COUPLERS 



51 



The matrix S may be written 

S = R X 5a + R b X Sb + R c X 5c + R d \ sd , (64) 

where the eigenvalues of S are 

^ = tti (65) 

(66) 

(67) 

--*+* (68) 

If the values of X fo , X f6 , X f „ , X fd are substituted in (65)-(68) and those 
are in turn substituted into (64) while the conditions indicated by 
(46) and (48) are imposed, the following result is obtained after some 
algebra 
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(70) 
(71) 

(72) 
(73) 

(74) 
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If, instead of the condition of (46), the one of (45) is imposed plus (48) 
then the following result is obtained 



S = 
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(75) 



where 



<Sl3 — 



S u — 



+ 



cosh y + 1 + sinh y + 1 cosh 7 I + sinh 7 Z ' 



cosh y + l + sinh 7 + Z cosh 7 Z + sinh 7 I ' 



where 



Z - 



v 



L u 



C + C M ' 



7 ™ 



7 = 



j(L u — Li 2 )co 

Z 



(76) 
(77) 

(78) 
(79) 
(80) 



The matrices of (69) and (75) with the aid of (70)-(74) and (76)-(80) 
give the frequency dependency of the coupling between the ports for 
the two types of directional couplers derived (matched at all frequencies) 



IV. EXPLANATION OF THE EFFECT 

Two sets of conditions have been derived for obtaining the directional 
coupler effect at all frequencies for transmission lines. Fig. 8 represents 
a pair of lossless lines c, d which are coupled from x = — I to x = 0. 
The coupled lines are connected to four (uncoupled) lossless transmission 
lines a, b, e, f, each of characteristic impedance R and each terminated 
in R. The matrix Z is the characteristic impedance matrix of the set 
of coupled lines (c and d in Fig. 8) and r its propagation matrix. 

Consider a case in which the eigenvalues of the matrix Z of the 
lines c, d satisfy 

R = Z = Z , 
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which implies 



Z = RI*, 



where I is the unit matrix. Assume that a pulse travelling down line 
a occupies at time t x the position shown in Fig. 8. Lines a and b together 
may be considered as a particular case of a set of multiple lines, in 
which both (Z ) ab and (r) o6 — the characteristic impedance and propaga- 
tion matrices of lines a, b — are diagonal because lines a and b are 
uncoupled. Thus it is convenient to think of the pulse traveling down 
line a as a vector of pulses traveling down the two lines a, b ; the second 
component of the vector which corresponds to line b being zero. When 
the vector of pulses arrives at the position x = — I, indicated by M 
in Fig. 8, the vector of pulses continue "seeing" the same characteristic 




UNCOUPLED 



Fig. 8 — Pair of coupled lossless lines c and d above a ground plane which has 
uncoupled terminated lines a, b, e, j of characteristic impedances R. This shows 
the progress of a pulse to explain the directional coupler effect in physical terms 
for a coupler with diagonal characteristic impedance matrix. 



impedance matrix RI as before and hence no reflection of the vector 
is created at M. 

In a vector formulation one speaks of reflections in a multidimen- 
sional sense. The voltage on line a may give rise to a reflected voltage 
on line a [self reflection] or to a reflected voltage on line b [mutual 
reflection]. If all the lines are self matched and mutually matched 
there will be no reflections whatever. Often a line might be self matched 
but not mutually matched; then no self reflection will occur, but a 
mutual reflection will. 

At time t 2 the second component of the vector pulse is no longer 
zero because the pulse on line c which is coupled to line d induces a 
pulse on line d as shown in Fig. 8. These component pulses will be 

* Although for simplicity in the explanations, parallel lines are assumed, this 
kind of coupler requires negative mutual inductances which are in general achieved 
with counter-wound helices. 
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distorted because continually varying portions of them travel at 
different speeds on the lines c and d* When they reach x = marked 
by N in Fig. 8, the vector again continues to see the same characteristic 
impedance matrix and hence no reflections are caused at N. Finally 
the pulses travel out on lines e and / and are dissipated at the resist- 
ances R. 

From this it is clear that the two ports at M are uncoupled and also 
the two ports at N are uncoupled. However a port at M is coupled to 
both ports at N and vice versa. All the ports are self matched. This 
explanation suggests a simple way of determining the conditions for a 
nonsymmetric coupler realized with transmission lines. Assume lines 
c and d are no longer identical but that the characteristic impedance 
matrix of the set is still diagonal 

\Z )u 

. (Z )i 

while the propagation matrix r is not. If the lines a and e have char- 
acteristic impedance (Z ) u and lines b and / have characteristic imped- 
ance (Z ) 22 and a, b, e, f are properly terminated, the coupled lines 
c, d will constitute a matched nonsymmetric directional coupler since 
the discussion above holds for this case without modification. 
A physical account of the directional coupler satisfying 

7 + = 7~ = 7, 



Z = 



«->P ; 



C + 2C M ' 



is as follows. 

Fig. 9 shows the same arrangement as Fig. 8. A pulse traveling 
to the right on line a is shown at t 4 t x . The second component of the 
incident pulse corresponding to line b is zero since lines a and b are 
uncoupled. As the vector of pulses reaches the position M, the vector 
of voltage is reflected according to a reflection matrix because for 
this case, the vector no longer sees the same characteristic impedance 
matrix in the transition from lines a, b to lines c, d. However, because 
line a matches line c there is no self reflection; only a mutual reflection 
appears on line 6. Besides the reflected pulse, an identical transmitted 
pulse appears on line d at time t 2 as indicated in Fig. 9. The appearance 

* The propagation matrix should not have equal eigenvalues, otherwise it will be 
diagonal which, together with a diagonal characteristic impedance, implies un- 
coupled lines. 
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b d f 

UNCOUPLED ] COUPLED COUPLED ] UNCOUPLED 

Fig. 9 — Same structure as in Fig. 8 showing the progress of a pulse for a 
coupler with scalar propagation matrix. 

of the mutual reflection indicates that the two ports at M are coupled. 

After t 2 , as the two forward pulses travel along lines c and d they do 
so at the same speed, without distortion and do not interact with each 
other since the propagation matrix is diagonal. As the two pulses 
arrive at point N they encounter a reflection matrix of opposite sign 
to the one they encountered at M. This means that the pulse on line c 
passes through N undisturbed but creates on lines d and / reflected 
and^transmitted pulses identical to the ones created at M but of op- 
posite sign. Likewise, the incident pulse on line d goes right through N 
(since each individual line is matched) creating transmitted and reflected 
pulses on lines e and c, but being cancelled on line / by the transmitted 
pulse created by the incident pulse on line c; thus, nothing comes out 
of line /. 

At time t a right after the reflection at N the situation is depicted 
in Fig. 9. After t a the reflected pulses are traveling to the left at the 
same speed undisturbed and undistorted on lines c and d. As they 
arrive from the right at point M the pulse on line d goes out line 6 
undisturbed but creating transmitted and reflected pulses on lines a 
and c. Likewise the pulse on line c goes out line a undisturbed creating 
transmitted and reflected pulses on lines b and d, but being cancelled 
on line a by the transmitted pulse created by the incident pulse on 
line d. This eliminates any delayed reflections on line a to the original 
incident pulse. The process continues in the same manner, the outgoing 
pulses on lines a and / always being such that they cancel. This means 
that the ports associated with lines a and / are uncoupled, but the 
ports of lines a, b and e are coupled. 

It is clear that what is necessary for directional coupler effect on this 
type of coupler is: all ports self matched, equal propagation velocities 
without attenuation or distortion. Hence it should be possible to 
realize a nonsymmetrical directional coupler of this type whose propaga- 



56 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 19G8 

tion matrix is a scalar matrix* by self matching its ports. This may 
be useful for interconnecting lines a and b of different characteristic 
impedances. 

Because the reasoning above was made in the time domain with 
pulses of arbitrary shape, the results hold for all frequencies. This 
is an example of the gain in insight owing to the vector-matrix for- 
mulation. 

V. EXTENSION OF THE THEORY 

5.1 Matrices jor a 2N-Port 

Let us generalize several concepts introduced in Section 2.1. Con- 
sider a (2N + l)-terminal network in which terminal 2N + 1 will 
be grounded and ports from terminals 1 through 2N to ground will 
be considered. Ports 1 to N will be considered input ports and Ports 
N + 1 to 2N output ports. Suppose the 2iV-port is characterized 
by A, B, C, D N X N matrices. (Extensions to 2iV-ports of the A, B, 
C, D parameters of a two-port). Assume the circuit is such that 

A - D, (81) 

A 2 - BC = I, (82) 

where I is the N X N unit matrix and A, B, C are N X N symmetric 
matrices which commute. 

By analogy with a multiple transmission line the characteristic 
impedance matrix Z and the propagation matrix r are denned so 
that they satisfy the following equations : 

A = cosh r, (83) 

B = sinh T, (84) 

C = Z" 1 sinh r. (85) 

Solving for r and Z 

r = cosh" 1 A = sinh -1 B, (86) 

Z = VBC 31 . (87) 

The N X N matrices r and Z will also be symmetric and commute. 
The matrix Z is the open circuit impedance matrix of that network 
which, when connected to the output ports N + 1 through 2N of 

* A scalar matrix is the unit matrix multiplied by a scalar. 
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the circuit whose A, B, C matrices are those of (86) and (87), will 
result in an open circuit impedance matrix of Z when the circuit 
is viewed at its ports 1 through N. This property is analogous to the 
one of the characteristic impedance matrices of a set of multiple 
coupled transmission lines. The matrix r, has virtually the same 
properties of the matrix rZ of a set of coupled lines (although the 
single quantity I looses its significance as a length in case the 22V-port 
is a lumped circuit). For instance, if n identical 22V-ports are cascaded 
the resulting 2iV-port has a propagation matrix equal to nT. 

The matrices r and Z may be expressed in terms of the 2iV X 2iV 
impedance matrix Z of the 2iV-port with the aid of (38) through (41) 

r = cosh -1 (Z n Z- 2 ;), (88) 



Z - Vg - Z 2 l2 ; (89) 

where the Z matrix is partitioned as follows : 



[ Zn Z 12 
Z 21 I Z 22 J 



(90) 



the submatrices Z n , Z 12 , Z 2I , Z 22 are N X N symmetric matrices 
and commute with each other. The characteristic impedance matrix 
may also be expressed in terms of the so-called open and short impedance 
matrices. If the iV-vector V\ and Ij denote the voltages and currents 
at the N input ports and V 2 and I 2 denote the voltages and currents 
at the output ports, then if I 2 = 0, that is, the terminals on the output 
ports are open then 

V, - AV 2 , (91) 

I, = CV 2 . (92) 

Solving for V 2 in (92) and substituting in (91) 

V, = AC-% . (93) 

which shows that the N X N impedance matrix Z 0c seen at the input 
ports is 

Z 0c = AC"'. (94) 

Now, if the output ports are shorted, that is V 2 = then 

V, = BI 2 , (95) 

I, = AI 2 , (96) 
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from which 

V x = BA-% . (97) 

That is the N X N impedance matrix with the output terminals shorted 
Z ae is 



Z 3c = BA" 1 . 

From Eqs. (94) and (98) it is seen that 



Z — V Z 5c Zoe . 



(98) 



(99) 



Equation (99) gives an experimental method of determining Z„ if 
it is known that the network satisfies equations (81) and (82), and 
the symmetry and commutativity conditions. 

5.2 Multi-port Circuits and Multiple Transmission 

It is often convenient to analyze some lumped or distributed (or 
combinations of lumped and distributed) systems as though they were 
multiple transmission lines using such concepts as reflection matrix 
and incident voltage. 

Consider the connection shown in Fig. 10. Each network is an (27V+1)- 
terminal network in which ports from each terminal to ground are 
made. Ports 1 through N and 1' through N' are considered input ports. 
Ports N + 1 through 2N and (N + 1)' through (2iV)' are considered 
output ports. The voltage vector at the junction B whose components 
are the voltages of nodes N + 1, N -j- 2, • • • , 2N to ground is denoted 
byV B . 

The vectors v + , v_ , i+ , i_ , called incident voltage, reflected voltage, 
incident current, and reflected current at the junction B (assuming 
the direction of propagation from left to right), are defined to satisfy 
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(2N)'i 



TV77777? 



Fig. 10 — Connection of two 2N ports in cascade used to define the reflection 
matrix. 
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V B = V + + V_ , (100) 

I* - U - L , (101) 

i + = Z"V + , (102) 

i- = Z -1 v_ ; (103) 

where I B is the AT-vector whose components are the currents flowing 
out of terminals N + 1, N + 2, • • • ,2N and into terminals 1', 2', • • • , N' ; 
Z is the N X N open-circuit impedance matrix of network M seen 
from the ports N + 1, N + 2, ••• ,2N with network Af' disconnected. 
The reflection matrix r^w is defined according to 

v + = r^-v. . (104) 

The matrix Y uu - satisfies the following relationship: 

IW = (Z'Z" 1 + ir'tZ'Z" 1 - I), (105) 

where Z' is the N X N open-circuit impedance matrix of network 
M' as seen from ports 1', 2', • • • , N' with network M disconnected. 
I is the N X N unit matrix. The indices MM' on IY^. indicate the 
direction of propagation from M to M '. If the indices are reversed 
the roles of Z' and Z are reversed, that is 

IV,, = (Z(Z')- 1 + I)(Z(Z')" 1 - I). (106) 

The transmission matrix T MM . is defined by 

V B = T AfJtf -v + . (107) 

Hence T MU - satisfies 

TW = I + r MM . . (108) 

5.3 Directional Coupler Equations 

Taking advantage of the derivations done for transmission line 
directional couplers and the analogies introduced in Sections 5.1 and 
5.2, it is possible to write without further work, the equations of a 
directional coupler having the same mathematical symmetry of a 
multiple transmission line directional coupler but which may have 
lumped components or combinations of lumped and distributed com- 
ponents. Suppose a four-port is characterized by its E matrix whose 
2X2 submatrices A, B, C, D have the form of the matrix K of equa- 
tion (1) of the Appendix and satisfy equations (81) and (82). Without 
any further work it can be stated that if the load impedance z and the 
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four port satisfy for some frequency the equation 

z = Z = Z , 



(109) 



the device will be a directional coupler for that frequency. 

This result is deduced from equation (45). Likewise from equation 
(46) it may be deduced that if the four port is such that 

(HO) 



7=7 



and 



z = Vz + z;, 



(in: 



then the device will also be a directional coupler for those frequncies 
for which (110) and (111) are satisfied. 

It is convenient at this point to exemplify with a simple lumped 
circuit. 

Consider the four-port lumped circuit shown in Fig. 11. The E 
matrix of the circuit of Fig. 11 may be calculated by cascading 3 sec- 
tions, the first and third containing only capacitors and the second 
containing the inductors and mutuals. By proceeding carefully much 
labor can be saved using the spectral sets given in the Appendix. The 
results are 
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Fig. 11 — Four-port lumped circuit that can be used as a directional coupler. 
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From (86) the hyperbolic sine and cosine of the propagation matrix is 



cosh r = - 



1 1 



[1 + S 2 C{L n + L 12 )] 



1 1 



+ 



1 -1 
-1 1 



[1 + S\C + 2CV)(L U - L 12 )], (113) 



sinh r 



=iri i 



i i 



S(L n + L i2 ) + - 



1 -1 
-1 1 



S(L U - L l2 ). (114) 



From (S7) the characteristic impedance matrix is found to be 



z„ = | 



1 1 

Li U 



(L„ + L 12 ) V 1 

[S 2 C(Ln + L 12 ) + 2]CJ "•" 2 



1 -1 
-1 1 



(, 



(Ln — L x2 ) 



v[S 2 (C + 2C„)(L n - L 12 ) + 2](C + 2CV) 
The condition expressed by (109) is, for this case, 
L n + L 12 



(115) 



z = 



[S 2 C(L XX + L 12 ) + 2]C 



L xx L n 



,[S 2 (C + 2CV)(L n - L 12 ) + 2](C + 2C U )> 
The second and third members of (116) imply 

[S\L 2 XX - L\ 2 ){C + 2C M ) + 2(L U + L 12 )](C + 2CV) 

= [S 2 C(L 2 n - L 2 2 ) + 2(L n - Lj^C. 
If this condition is to be satisfied at all frequencies then 
(L 2 , - L\ 2 ){C + 2C M ) 2 = (L 2 1X - L\ 2 )C 2 , 
2(L„ + L X2 )(C + 2C M ) = 2(L n - L X2 )C. 



(116) 



(117) 

(118) 
(119) 



Both (118) and (119) are satisfied for the following choice: L u = —L 12 , 
C = 0. Thus the circuit of Fig. 11 with C = and L xx = —L i2 (per- 
fectly coupled counterwound inductors) is a directional coupler at all 
frequencies, provided it is loaded at all ports with the impedance 



z = 



L xl 



.4L„C*S J + 2CJ * 



(120) 
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The impedance z is frequency dependent. The voltage ratios will be 
given by (76) and (77),* that is 



Sl * ~ 2 I 1 + 4C U L U & + 1/ ' 



(121) 



S -=\\ l - ,C M L^ + & (122) 

5.4 Equations in Terms of A, B, C, D Matrices 

It is often convenient to express the equations of a directional coupler 
in terms of the A, B, C, D matrices directly instead of the Z and r 
matrices. For this purpose assume a lossless reciprocal four-port is 
characterized in terms of its A, B, C, D matrices which are of the form 
of the matrix K of equation (150) of the Appendix. Assume A = D. 
The condition A 2 — BC = I is automatically satisfied if the circuit 
is reciprocal. A, B, C, D commute, since they have the same eigen- 
vectors. Because all matrices commute they may be treated without 
ambiguity as scalars. The open circuit impedance matrix is 



C = 



AC 



AC 



•} 



(123) 



Suppose the ports are loaded with equal impedance z. The impedance 
matrix < normalized with respect to the matrix zl is 

The eigenvalues of < n are 

A„ = ^T" 1 , (125) 

A» - A * C T 1 , (126) 

K = ££± , (127) 

a« - ^^ ; (128) 

where A + , A~, C* ' , C~ are the eigenvalues of A and C associated with 
the sum and difference modes. (See Appendix.) The reflection matrix 

* Because for lumped elements r corresponds to rl in using the formulas derived 
for distributed elements for circuits with lumped elements one should take |»1 
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(or scattering matrix) is 

r K - «„ - i)« B + 1)- 1 , (129) 

which may be written 

r «= R -m + R '^TT + R 't!TT + ^m ; (130) 

where R a , R 6 , R c , R d are the members of the spectral set of {„ and 
which are given by equations (165) to (168) of the appendix. 

To find the condition for self match at port 1 (which because of the 
symmetry gives the condition of self match at any port) the eigenvalues 
given by equations (125) to (128) are substituted in (130) and the 
upper left corner of r R is equated to zero. After some algebra this yields 

[(a- + zcy - i][(A + ) 2 - (zcy - i] 

+ [(A + + zCy - 1][(A-) 2 - {zC~f - 1] = 0. (131) 

Equation (131) is a quartic in z which may be rewritten 

z\C + C~) + z\C + A~ + C~A + ) 

- z(A + B~ + B + A~) - B + B~ = (132) 

The solutions of Equation (US) give the values of the impedances which 
will match the four ports in terms of the eigenvalues of the matrices 
A, B, C. Although a quartic algebraic equation can be solved in terms 
of the coefficients, the solution is extremely cumbersome algebraically 
and it would be very difficult to see the effect of varying the quantities 
A + , A~, B + , B~, C + , C~. A sounder approach is probably to look at 
particular simple cases. For instance if 

C + = and B~ = 0, (133) 

Equation (132) reduces to 

z\C'A + ) - z(B + A~) = 0, 
whose solutions are 



z = and z = ^/§=jr- (134) 

A second possibility is 

C* =0, A + = 0. (135) 
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Eq. (132) then reduces to 

-z{B+A~) - B + B~ = 0, 
whose solution is 

z = -~ (136) 

Obviously there are many possibilities. Some of the solutions are not 
immediately apparent. For instance from equations (46) and (48) if 

A + = A", (137) 

then 



z = 4 - S /p3 (138) 



'('- 



is a root of equation (132). This fact can only be seen after a good deal 
of algebra, for this reason it is convenient to express (132) in different 
ways so that different possibilities may be "seen." With this in mind 
equation (129) may be written 

r R =1- 2(C+ir 1 . (139) 

Using the spectral set of <„ the following alternative expression for 
the condition for the self-match of all ports is obtained 

12 12 



4 A+ + 1 "4 A+-1 

zC + + 1 zC + + 



, I ? 1 2 

+ 4 A- + 1 +4A--1 ' 

zC~ + zC~ r 

which after some algebra may be written 

( A + + zC + )zC + , (A~ + zC~)zC- _ (im 

(I + + zcy - i + (a- + zcy - i " *■ U4U; 

It is simpler (although not trivial) to verify that the conclusions asso- 
ciated with equations (137) and (138) are true from (140) than from 
(132). 

The reflected voltages caused by an incident voltage at port 1 may 
be obtained from equation (130). 
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- 1 


1] 2[(A- + zcy - 


- 1] 


i * c " 




1 (A - + zC") 2 - 1 ' 




zC~ 





(141) 
(142) 
(143) 



Equations (123) through (143) all are good for any four-port, whether 
lumped, distributed, or made with combinations of lumped and dis- 
tributed elements. 

5.5 A Degenerate Situation: 

Consider the circuit of Fig. 12. Notice the structure is not physically 
symmetrical. The A, B, C, D matrices of the circuit are 



|_c I dJ r= Ly I iJLo 1 1 _ 



YZ +I_ 



(144) 



where 



Z - * l 2SL„ , 



Y = 



-1 
1 



28C, 



Although in equation (144) A and D are apparently not equal, it turns 
out that Y and Z are orthogonal and therefore YZ = 0. Thus (144) reads 



A 



B 
D 



(145) 



The matrices A, B, C, D commute and satisfy equations (81) and 
(82). Thus, although the structure is not physically symmetrical, it 
is electrically symmetrical. When one attempts to use equation (87) 
to determine Z one finds that the matrix C" 1 does not exist because 
C = Y is singular. Since Z does not exist, equations (48), (70), and 
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Fig. 12 — Lumped directional coupler which, when connected to a resistive load, 
exhibits directional coupler effect of all frequencies. 
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(71) cannot be used. Equation (132) can be used instead. The eigen- 
values of A, B, C are : 







A + 


= 1, 




A" 


= 


1 
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= 4SL 


ii > 


B~ 


= 











C + 


= 0, 
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4SCV 




Thus (132) 


reads 
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(146) 

The load impedances are frequency invariant, which indicates that 
the circuit may be matched with a constant at all frequencies and 
should exhibit directional coupler effect at all frequencies when loaded 
with a positive resistance of value \/L u /C ' u . Using equations (141)— 
(143) the voltage ratios are found to be 

V2- _ 2S'yL n C M (147) 

Vi+ 1 + 2SVL U C U ' 

¥*=■ _ 1 (148) 

V 1+ l + 2SVL U C M ' 

= 0. • (149) 



v 1+ 



Equation (149) corroborates that the coupler exhibits directional coupler 
effect at all frequencies. This example illustrates the use of the direc- 
tional coupler equations in terms of the A, B, C, D matrices. 

5.6 Lumped and Distributed Elements 

The formulation that has been developed allows the handling of 
circuits with both lumped and distributed elements without any changes 
because the formulas are good for "black boxes." For example, for 
the circuit shown in Fig. 13, the total E matrix is found by multiplying 
the individual E matrices of the sections. The E matrix of section P 
or T is given by equation (145) while that of Q is given by equation 
(29). Once the total E matrix is known, it is partitioned into A, B, C, D 
matrices and equation (132) applied to determine the proper z for 
terminating the coupler. When the coupler is thus terminated, equa- 
tions (141)-(143) yield the voltage ratios. 
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Fig. 13 — Directional coupler configuration containing both lumped and dis- 
tributed elements. 

In general, the algebra will get quite unmanageable and one will 
have to resort to numerical calculations on a digital computer at 
discrete frequencies. 13 The z may be found by numerically solving 
the quartic equation (132) at a set of discrete frequencies and then 
realizing it as a driving point impedance through successive approx- 
imations, or some similar procedure. The processes of normalizing the 
impedance z and of making frequency transformations can be used 
very effectively in the realization of directional couplers of this sort. 



VI. CONCLUSIONS 

Although strictly speaking all physical devices are distributed in 
space and thus, in general, have transcendental transfer functions for 
certain frequency regions, it might be possible to model the devices 
accurately enough with conventional ideal lumped elements, or more 
generally with elements having given frequency curves, which may 
be given analytically or numerically. In this paper we give a matrix 
theory for lumped and distributed circuits, keeping this fact in mind. 

By using matrix formulation and treating the circuits as black boxes, 
it is possible to extend the classic theory of stripline directional couplers 
to more general circuits while still keeping many of the concepts (such 
as even and odd characteristic impedance) that have been found useful. 

The paper makes evident the fact that the concepts of odd and even 
mode arise because of the special symmetry of the matrices and that 
they correspond to their eigenvectors and eigenvalues. We indicate 
in the appendix that when such special symmetry is lost, the odd and 
even modes are also lost and it might be necessary to introduce one 
set of modes for the currents and another for the voltages. This fact 
is not simple to see without the matrix formulation. 

Bj r thinking in vector-matrix terms, we explain the directional coupler 
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effect and gain considerable insight which is useful for realizing direc- 
tional couplers lacking double symmetry. 

We have given general equations for analyzing and designing black- 
box directional couplers in terms of the characteristic impedance and 
propagation matrices, and in terms of the transmission A, B, C, D 
matrices. The latter may be necessary to analyze circuits whose char- 
acteristic impedance matrix (hence even or odd characteristic imped- 
ances) does not exist but which have a scattering matrix. 

The topic of the actual design of directional couplers with lumped 
or with lumped and distributed elements, and more specifically the 
design of multisection directional couplers using computer aids, is not 
treated because it is the subject of a forthcoming paper. 



APPENDIX 



Here are some spectral properties of the principal matrices in the 
paper for ease of reference. 14 
The symmetric matrix 



K = 



_k 12 fen.. 



has eigenvalues 



Ai — A'n T~ "-is 

A2 = K\\ /Cis 



and normalized eigenvectors 



U x = 



1 



V^ 



u, = 



V~2 



1], 
-U 



This can be seen by verifying the following identity 



1 



1 1 
l1 -1. 



\/2 
The spectral set of K is: 



fc„ + k i2 

Ku fci2_ 



1 1 
1 -1. 
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k\2 ">ii_ 
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R2 = TJ2U2 j R2 = 
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(150) 

(151) 

(152) 

(153) 

(154) 

(155) 
(156) 



DIRECTIONAL COUPLERS 



69 



Hence an analytic function / of the matrix K may be written 

/(K) = R,/(X,) + R 2 /(X 2 ) (157) 

For example 

cosh *» kl2 

._/C 12 ACii_. 

= r|[cosh(fc 11 +fc 12 ) + cosh(fc 11 -fc 12 )] £[cosh(fc n +fc 12 )-cosh(fc n -fc 12 )]~| 

_^[cosh(/.- 11 +fc 12 )-cosh(fc 11 -A- 12 )] ^[coshCfcu+fcj^+coshCfcn— fc 12 )]j 

The quantities associated with the vector U l are called the "even mode" 
or "sum mode." The quantities associated with the vector U 2 are 
called the "odd mode" or "difference mode." 

In our main paper, the eigenvalues of the matrices Z and r which 
are in the form of equation (1), are denoted by Z\ and 7 + for the sum 
mode and Z~ and y~ for the difference mode. 

The partitioned matrix 



M = 



M„ 



_M 1S 



M i: 



M„_ ' 



(158) 



where M n and M 12 are the following 2X2 symmetric matrices 







M„ = 


J o 
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7 
5 


5 

7J 


> 




(159) 


has the following eigenvalues 








X„ = (a + /3) + (7 + 5), 






(160) 


X, - (a + /S) - (y + 8), 






(161) 


K = (« - 0) + (y - 8), 






(162) 


\ d = (a - 0) - (y - 8), 






(163) 


and the corresponding eigenvectors 
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• (164) 
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This can be verified by a matrix multiplication similar to that of equa- 
tion (154). 
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The spectral set of M is 
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(165) 



(166) 



(167) 



(168) 



Any analytic function / of M may be written 

/(M) = R„/(X a ) 4- R„/(X 6 ) + RJ(X C ) + R„/(X„). (169) 

Concerning the so-called "modes of propagation" of a set of two 
coupled lines, when the r matrix of the two lines has the double sym- 
metry exhibited by the matrix K of equation (1), the eigenvectors 
of the matrix are those given by equation (4). Since the matrix is 
symmetrical, the eigenvectors of the transposed matrix r' are the 
same; therefore, one may speak of the "sum mode voltages and cur- 
rents" and "difference mode voltages and currents." However, it might 
happen that r does not have the form of K in equation (1). Then the 
concept of sum and difference modes disappear because the eigenvectors 
that will result will not be quite so simple. If he wishes, one may then 
speak of "first mode" and "second mode," associating each mode with 
each eigenvector. 
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If the matrix r is symmetrical the voltage and current modes will 
coincide. However, if r is not symmetrical then the eigenvectors of 
r' will not be the same as those of r. It will be necessary to speak of 
"first voltage mode," "second voltage mode," "first current mode," 
and "second current mode" because the voltage modes will differ from 
the current modes if r is not symmetrical. The eigenvalues of a matrix 
and its transpose are always the same, hence no distinction is necessary 
for the propagation constants of the voltage and current modes. 
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